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Abstract — The present study is concerned with water wave scattering by a porous bottom of varying depth in presence of a 
discontinuity at the upper surface of the ocean. Assuming linear theory and using a perturbation technique in 
conjunction with the Green's integral theorem, the first order corrections to reflection and transmission coefficients are 
obtained. By considering two different shape functions describing the bottom undulation, the effect of porosity is 
investigated numerically and these results are depicted in a number of figures against the wave number of the incident 
wave. 
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I. INTRODUCTION 

The problem of water wave scattering due to uneven bottom is an interesting topic of research over a few decades. 
A class of water wave scattering problems in presence of a small obstacle situated at the bottom was considered earlier by 
many researchers such as Lamb [1], Kreisel [2], Davies [3-4]. Mei [5] considered the problem of reflection of water waves 
by periodic sand bars and explained a theory on strong reflection i.e Bragg reflection induced by the sandbars themselves. 
Mandal and Basu [6] extended this problem to include the effect of surface tension at the free surface for an obliquely 
incident surface wave train and by employing a perturbational analysis directly to the governing problem to obtain reflection 
and transmission coefficients up to first-order in terms of integrals involving the shape function describing the bottom 
topography. Evans and Linton [7] considered the problem of water wave scattering by uneven bottom using step 
approximation and introduced scattering matrix. Ultimately they reduce this problem to a problem of scattering by a 
discontinuity at the upper surface in a uniform finite depth water and employed residue calculus technique of complex 
variable theory to obtain the reflection and transmission coefficients of the incident wave. Martha and Bora [8] also 
investigated a scattering problem in presence of a small bottom undulation and by using perturbation technique they obtained 
the reflection and transmission coefficients upto the first order. Recently Mandal and De [9] investigated a wave scattering 
problem by a free surface discontinuity together with a small bottom undulation. Using Green's integral theorem they 
obtained reflection and transmission coefficients up to first order in terms of computable integral. 

Problem of water wave interaction with a rapidly varying porous bottom is another important aspect in this 
context. The flow of fluid into the porous media leads to different phenomenon like wave energy dissipation, damping etc. 
Water wave interaction with the porous media was studied earlier by many scientists like Chwang [10], Chakrabarti [11], 
Jeng [12] and many others. Zhu [13] studied the wave propagation problem within porous media on an undulating bed by 
employing a Galerkin eigenfunction expansion technique and investigated the wave reflection coefficient numerically. Gu 
and Wang [14] also investigated the problem of water wave interaction with rigid porous seabed both theoretically and 
experimentally. Mase and Takeba [15] focused on the Bragg scattering of gravity wave over a porous sea bed. Also the work 
of Silva [16] discussed on water-wave reflection and transmission problem where a porous medium was assumed to lie on a 
rapidly varying sea-bed. Recently, Martha and Bora [17] considered the problem of oblique water wave scatting by a small 
patch of sinusoidal bottom undulation over a rigid porous bottom and obtained the reflection and transmission coefficients 
by employing fourier transformation technique. 

In the present article, we consider a water wave scattering problem by a rigid porous bottom of variable depth in 
presence discontinuity at the upper surface of the ocean. The problem is mathematically formulated in case of normally 
incident wave. The discontinuity at the upper surface may be presented in different form such as two different vast floating 
inertial surfaces of non-interacting materials or large broken ice sheets in the Antarctic regions etc. Here the upper surface of 
the ocean is assumed to be covered by two vast floating inertial surfaces of different thickness which produces a 
discontinuity at the upper surface of water. The effect of porous bottom is taken into account by assuming that the ocean 
bottom is composed of some specific type of rigid porous material with small undulation which characterized by a known 
shape function. The bottom boundary condition in this case usually involves a known porosity parameter G representing the 
typical of porous bottom. The wave train is incident from negative infinity direction which is partially reflected and partially 
transmitted along the surface of the ocean. The presence of small undulation at the bottom suggests that a perturbation 
technique can directly be applied to the governing boundary value problem and the associated boundary conditions. By 
applying perturbation analysis in terms of a small parameter e characterizing the bottom undulation upto order one, two 
different boundary value problems(BVP) are generated. They are renamed as BVP-1 and BVP-2 respectively. The BVP for 
the zero order potential function represents a scattering problem in presence of a discontinuity at the upper surface over a 
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uniform porous bottom. This problem is solved here by using residue calculus technique followed by Evans and Linton [7] 
and the zero order corrections to the reflection and transmission coefficients are obtained for the purpose. Now without 
solving the BVP-2 for the first order potential function, the first order corrections to the reflection and transmission 
coefficient are obtained by using Green's integral theorem. The analytical expressions for these coefficients are presented in 
terms of computable integrals involving the bottom shape function, zero order potential function and an additional term 
containing G. We choose two different bottom profile functions, one is exponential and other sinusoidal to describe the 
undulations at the ocean bottom. These type of bottom profiles can be observed in different regions of an ocean due to 
growth of sand ripples, sedimentation of different porous materials at the bottom and other naturally occurring phenomenon. 
For these type of bottom profiles, the first order corrections for reflection and transmission coefficients are depicted 
graphically against the wave number of the incident wave in a number of figures. The effect of porosity on the first order 
reflection and transmission coefficients is investigated in presence of surface discontinuity by taking different values of the 
dimensionless porosity parameter in different cases. 

II. MATHEMATICAL FORMULATION OF THE PROBLEM 

A right handed regular cartesian co-ordinate system {x, y) is chosen in which y — represents the undisturbed 
upper surface of water and the y-axis is taken vertically downward into the fluid region. The upper surface is covered by two 
floating inertial surfaces of uniform area densities Eip and £ 2 p (p being the density of water). One of the inertial surface 
occupies the region y = 0,x < and the other occupies the region y — 0,x > 0. The wave field is incident from x -> -co 
and the direction of the positive x-axis being opposite to the direction of the incoming incident wave field. Assuming linear 
theory and irrotational motion, the velocity potential describing the time-harmonic motion of angular frequency a, can be 
represented by Re[cp(x,y)e~"' t ], where <p satisfies the following two dimensional Laplace equation: 

—■ + d -^r = 0, onO < y < h + ec(x), -oo < X < 00. (1) 

The upper surface boundary conditions are prescribed by 

K l( p + d J-= ony = 0,x<0, (2) 

K 2( p + ^- = ony = 0,x>0. (3) 

This produces a discontinuity in the upper surface boundary conditions at the point (0,0), where K 1 — 
,K 2 — , £1,62 < -j and K — — (g is the acceleration due to gravity). 

1— £\K 1—£2K o g 

The bottom boundary condition in case of uneven porous bottom (Martha, Bora, Chakrabarti [17])is given by 

d -^--Gcp = 0ony = h + ec(x). (4) 

Here y — h + ec{x) describes the bottom of the ocean with variable depth. The positive number £ signifies smallness of the 
undulation at the bottom and y\ is the outward normal to the ocean bed. The bottom undulation is characterized by a known 
shape function c{x) with the property that c{x) -» as x -> +co ensuring that far away from the undulation the ocean is of 
uniform finite depth of h below the mean free surface. 
The far field behavior of <p(x, y) is described by 

0(x , y) J^7o 2 W as,- cc, 

w y {(e lk » x + Re- lk « x )ipl(y) asx - -oo, 
where 

*Po(y) = N^(coshk (h - y) - - sinhfe {h - y)), 

"a 

^o(y) = W 2 (coshs (/i - y) - f sinhs (h - y)) 

3 

2*| 

2/c (G-G 2 h+/c^)-2G/c cosh 2k h+(ki+G 2 )sinh 2k h 

3 

' V - I =■ 

I2sq(G— G 2 h+SQh)— 2Gs cosh 2s /i+(sg +G 2 )sinh 2s /i 

Here, e lk " x xp} ) (y) represents the incident wavefield, R and T are respectively the unknown reflection and transmission 
coefficients(complex) to be determined. k and s are the unique breal positive roots (Cf. Mclver [18]) of the following two 
transcendental equations in terms of X: 
(1 + ^)tanhlh = #! + G, 

A 

(1 + ^)tanhAh = K 2 + G. 



and JVn 
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As £ is small, the approximate form of the bottom boundary condition (4) after neglecting 0(e 2 ) terms is 



III. METHOD OF SOLUTION 

As e is small, the app 
given by 

2i-G*- e (i{cM2^} + Gc(x)^) + O( £ 2) = ony = h. (6) 

The form of the approximate boundary condition (6) and smallness of the parameter £ suggests that (p, R and T has the 
following perturbational expansion in terms of £: 

<p(x,y;e) = <p (x,y) + £<pi(x,y) + 0(£ z ), (7) 

R(e) =R + efii + 0( £ 2 ), (8) 

7(£) = T + e7i + 0(£ 2 ). (9) 

Substituting (7)-(9) in (l)-(5), we find that (p and <p 1 satisfies two different boundary value problems. 

BVP-1 

The function cp (x, y) satisfies the following boundary value problem: 

?!^> + £f% = in < y < ft, -oo <x<<v, 

K 1 cp +Y^- = Oony = 0,x<0, 

K 2 cp +Y 1 = 0ony = 0,x>0, 

^ - Gcp = ony = h, 

a, ( x v ) _» i T o eis ° Xj Po(y) as x -» =o, 

Here i? and T stands for the reflection and transmission coefficient of zero order respectively. The above BVP of (p (x, y) 
describes a scattering problem by a discontinuity at the free surface boundary condition over a uniform porous bottom. 
Following Evans and Linton [7], using residue calculus method of complex variable theory the zero order reflection and 
transmission coefficients are respectively given by 

R = k J>ZE! e 2ia (10) 

T = .i*°- e '(«+« (11) 



-1^ _hn-lA 



Where, 

B = 5S?-i[tan- 1 (^-)-tan- 1 (^)], 
j8 = S?-i[tan- 1 (f2.)-tan- 1 (2.)]. 

Here, k n and s n (n — 1,2,3, ) are given by positive real roots of the following two transcendental equations in terms of 

A: 

(A - ^)tanAft + K t + G = 0, 

(1 - ^Jtanlft + K 2 + G = 0. 
BVP-2 

The first order potential function (p 1 (x, y) satisfies the following BVP: 
^-4± + 2-£l =0 in < y < ft, -oo < x < oo, 

K l<t>l +Y 1 =Oony = 0,x < 0, 

#201 + J 1 - ony = 0, x > 0, 
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2*1 - C0! = -{c(*)^} + Gc(x) d -^ony = h, 

dy ri dx l v ' dx ' v ' dy * 

, , . f7ie ,s °*i/ , o(y) as * ""> °°' 
I^Rje ° Vo(y) asx -> — qo. 

The above BVP for 0! represents a radiation problem in presence of an upper surface discontinuity over an uneven porous 
bottom.i?! and 7\ are the first order reflection and transmission coefficient respectively. It may be noted that bottom 
boundary condition for (pi involves the porosity parameter G and the known bottom shape function c(x). This ensures that 
the bottom undulation over porous bottom has an effect on the reflection and transmission coefficients at first order. 

The eigenfunction expansion of the zero order potential function cp (x,y) is given by 



T e is ° x xl>l(y) + S?-i B n e- S " x xti(y) as x > 0, 

(e"" + R e- ik » x )^(y) + E?-i A n e Kx Vn(y) as x < 0, 






where v4 n and B n (n = 1,2, ) are the unknown constants, ipniy), 0„(y)(n = 1,2, ) are the orthogonal depth 

eigenfunctions given by 

ipniy) - fy}(cosk n (7i - y) - — sinfc n (7i - y)), 

0«GO = N n 2 (coss n (h-y)-f sins n (h-y))(n = 1,2 ) 



where, 



n! 



*J2k n (-G+G 2 h+k 2 h)-2Gk n cos2k n h+(k 2 -G 2 )sm2k n h' 

3 " ' 

N n 2 = . , , 2s " , , (n = 1,2,3 ). 

V2s„(-G+G 2 h+s|h)-2Gs„cos2s„h+(sJ- G 2 )sin2s n ft 

Now using the matching condition of <p (x,y) at x — and orthogonality of the depth eigenfunctions produces the 
following system of linear equations 

y w B n NJ;coss n h _ r W 2 coshs h 

s n~k m ls 0^~k m 

(13) 



^,»>M j = N i coshkh f^o i_l (m = ! 2 } 



The unknown constants A n , B n (n — 1,2,3, ) are to be eliminated numerically from (13) after truncating the infinite 

sums upto the desired accuracy. 

IV. REFLECTION AND TRANSMISSION COEFFICIENT OF FIRST ORDER 

The first order corrections for the reflection and transmission coefficients can be obtained by using Green's 
integral theorem for the two potential function <po(x, y) and 0i(x, y) as 

& (0o0iv " <M>ov)<« = 0, (14) 

where v is the outward normal to the line element dl and the contour L is specified by the lines 

y = and - X < x < X, 

x = ±X and < y < h, 

y -h and - X < x < X{X > 0). 

The outgoing nature of the potential function (p (x,y) as x -> oo and the boundary conditions for <po(x, y) and <p\(x, y) 
ensures that only contribution to (14) comes from the bottom part y — hand — X < x < X and x = —Xand < y < h. 

Therefore, we obtain, 

2ifc fii = C(c(*X0iL(x,ft) - G 2 <pl(x,h))dx. (15) 

The first order transmission coefficient can be obtained by considering 0o(— x, y) and cp 1 (x,y) in (14) and arguing in the 
same manner as above, we obtain, 
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2is 7\ = - 5* m {c{x){<p 0x {x,h)<p 0x {-x,h) + G 2 <p {x,h)<p {-x,h))dx. (16) 

The integral representations of R 1 and 7\ given by (15) and (16) contains an additional term of porosity parameter G. 
However, in case of non-porous bottom(G=0), these forms are similar to the results obtained by Mandal and De [9] for the 
same type of problem. 

V. NUMERICAL RESULTS 

In this section, two different bottom shape functions are considered and numerical results for the first order 
reflection and transmission coefficients are explained briefly. 

Case-1: The first bottom profile function is taken as 

c(x) — Coe - * 1 '*' — oo < x < oo. 

The above c(x) represents exponentially decaying topography over porous bottom. The expressions of Ri and 7\ 
for this c(x) are given in the Appendix. The first order reflection and transmission coefficients are computed numerically for 
different values of Kh by taking c /h = 0.1, flh = 1.0, £±/h — 0.01, s 2 /h = 0.02 and the values of the dimensionless 
porosity parameter is taken as Gh = 0.00,0.05,0.1 respectively. These numerical results for Ri and 7^ are depicted against 
Kh in figure- 1 and figure-2 respectively against the wave number of the incoming wave. In the figure- 1 the amplitude of Ri 
is plotted against Kh for different values of the porosity parameter Gh. It is interesting to note that for each value of Gh, the 
graph of | /J-l | first increases with Kh up to a maximum value and then gradually decreases i.e a Bragg resonance is observed 
here in all cases. 
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Fig-1 

This may be attributed due to exponentially decaying bottom undulation and the nature of the incident wave field. 
These values of \R^ are observed to be rapidly decreasing with Kh which is very prominent for Kh > 1.256 because far 
away from the bottom undulation the ocean is of uniform finite depth h below the mean free surface. Moreover, for each Gh, 
occurrence of zeros of IRj Ifor certain values of Kh implies that the sinusoidal bottom does not affect the incident waves at 
first-order for certain frequencies. 
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Fig. 3 

Another interesting feature of figure- 1 is that as numerical value of the porosity parameter Gh increases, the 
highest amplitude of R 1 also increases with Gh. This can be explained as ^ in (15) contains c(x) and an additional term 
containing G 2 . Moreover the ocean bed is assumed to be composed of specific type of porous material with porosity 
parameter G along with the free surface discontinuity has an effect on the reflection and transmission coefficient at first 
order. 

In the figure-2 the amplitude of 7\ is plotted against Kh for different values of Gh. For each Gh, |7\| picks up a 
highest value for a certain Kh and gradually decreases with Kh which is prominent for Kh > 1.5898 in all cases. The 
highest amplitude of 7\ is found to be increasing with the increasing value of Gh. Also from the figure- 1 and figure-2, it is 
observed that figures of \R^\ and |7\| agrees with the results obtained by Mandal and De [9] in case of non-porous 
bottom(G/i = 0). 
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In figure-3, the amplitude of R 1 is shown against Kh for different values of jxh when c /h = 0.1, s^/h = 
0.01, s 2 /h = 0.02 and value of porosity parameter is taken as Gh = 0.05. It is seen that for each value of [ih, the nature of 
the graph of |/?i | is almost similar to case of figure- 1. In this case, the maximum amplitude of ^ is found to be decreasing 
with the increasing value of pJn.. This fact can be explained according to the form of R r given in the Appendix which shows 
that R 1 is inversely related to fih. In figure-4 the same figure is plotted in case of non-porous bottom (Gh = 0) and this 
matches with the figure obtained by Mandal and De [9] for the exponential bottom profile. 
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Case-2: The sinusoidally varying bottom profile is taken as 



c(x) = c sinlx — — < x < — — 

A A 



— otherwise , 
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where ripple amplitude c and ripple wave number X. The expression of R 1 is given in the Appendix for the above c(x). By 
(17) we get 7\ = 0. In figure-5, \R X \ is plotted against Kh by taking Gh = 0.00,0.05,0.1 when ejh = 0.01,£ 2 //i = 0.02 
and Xh — 1. This figure shows the changes of \R 1 \ for different value of the porosity Gh. The graph of 1^1 is found to be 
oscillatory in nature for each Gh and the oscillation gradually decreases with Kh which is prominent for Kh > 1.1. The 
oscillation of \R 1 \ picks a maximum value of for certain Kh hence a Bragg resonance is observed here in all cases. It is also 
found that as the value of Gh is increasing, the highest amplitude of R 1 is also found to be increasing and has shifted to the 
left. This may be attributed due to a sinusoidally varying porous bottom and the presence of surface discontinuity at the 
upper free surface. In figure-5, for case of non-porous bottom(G/i = 0), the graph of R 1 is almost similar to the results 
obtained by Mandal and De [9]. 
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Fig-6 

The figure-6 is plotted for \Ri | against Kh by taking two different values of the ripple amplitude c /h — 0.1,0.2 over porous 
bottom(Gh=0.05) while £i/h = 0.01, E 2 /h = 0.02 and Xh — 1. It is observed that overall values of \R^\ are increasing with 
the value of c /h. The same figure for \Ri\ in case of non-porous bottom(G/i = 0) is plotted in figure-7, agrees with the 
result obtained by Mandal and De [9] . 
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VI. CONCLUSION 

A simplified perturbation technique and Green's integral theorem are employed to obtain the first order reflection 
and transmission coefficient. From the numerical results and figures, it is observed that for a sinusoidal patch and 
exponentially decaying porous bottom with an upper surface discontinuity has a significant impact on the reflection and 
transmission coefficients at first order in different cases. The analytical and numerical results agrees with the results obtained 
by Mandal and De [9] in the case of non-porous bottom(G = 0). 

Problems of this type relating to porous bo ttom is found to be important in different branch of coastal science and 
engineering over a few decades. Understanding the effect of porous bottom on the wave characteristics in presence of a 
surface discontinuity over uneven bottom has gained considerable amount of interest to the scientists due to its various 
potential applications in different research areas of marine science and oceanography. 
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APPENDIX 

Case-1 (for sinusoidally varying bottom topography c{x) = c sinlx): 



R ^2k lHN ^ mo+G ^w^{^ me i - 1} ) + #-L ( ( - 1)me A - i ) +2R " (fe o 



((-l) m - 1) 1 + 

mn fX 

k„A„Nil I mn(k n +ik a ) \ \-l« k„A„Nil I mnjk n -ik ) 



+ G 2 ) ^^_ 1_ }+ ( \ k n A n e k "*Nl) 2 -G 2 Q k n A n e k "XN*) 2 )sirdxdx 



Z- hl=1 A 2 + (k n + ik yv ' J ° °Z- ln=1 A 2 + (k n -ik yv ' J 

mn: 

+ J ((2, ^Bne'^N 2 ) 2 -^ i B n e- s n"N 2 ) 2 j smAxdx - ((s 2 + G 2 )T 2 ^ _ — 2 ((-l) m e~^ - 1) j 
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-2^ro ^ n= a z + (5n _ iSQ y {(-ITe-^- - l) - 26'ToNJ ^ p - g- + ^ ((-l)-e^ - l)] 



and 

7\ =0. 

Case-2 (for exponential bottom topography c(x) = c e~^ x ^): 



H+(k n +ik a ) fi+(k n -ik a ) 

-2G 2 < S? =1 -^g- - ZG 2 ^ 1 ^ E^i u+ f; W "., : + / " ((E?=i s n B n e-^JV n 2 ) 2 

-G 2 (E^i B n e-^Nif)e~^Ax + (iV 2 ) 2 T 2 f^J - 2is N§T Tn=i '"*""; : 



and 



-2G 2 T N$YIn=i , J 



+G 2 T N 1 N 2 { .} , + .f° . } 

j_ *ti |\r2 Vw ™nkn"n -j rtI vioo H n k n N n _ y*oo Hnkn"n 

— ISnl niVn 2j1 " ; iKnlVn /.i — r fin /.i — 

G 2 T N§ Y.i . ^"5 + GX 1 Si ./""" + i? Ei ., g """ 

u u l is -k n -ii n ' J1 i/c +s„+^ u ' J1 ik -s n -ii 

-So ffi?-i fcn4.e-<*« + '0*iV B 1 )(E?-i k n B n e- s n*Nt)dx] 
+G 2 /; (S^i ^e-fc^^JV^CS^i S n e- S »^ n 2 )dx]. 
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